Phases of a 2D Bose Gas in an Optical Lattice 
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Ultra-cold atoms in optical lattices realize simple, fundamental models in condensed matter 
physics. Our ®^Rb Bose-Einstein condensate is confined in a harmonic trapping potential to which 
we add an optical lattice potential. Here we realize the 2D Bose-Hubbard Hamiltonian and focus 
on the effects of the harmonic trap, not present in bulk condensed matter systems. By measuring 
condensate fraction we identify the transition from superfluid to Mott insulator as a function of 
atom density and lattice depth. Our results are in excellent agreement with the quantum Monte 
Carlo universal state diagram, suitable for trapped systems, introduced by Rigol et al. (Phys. Rev. 
A 79, 053605 (2009)). 



The Bose-Hubbard (BH) Hamiltonian realized by 
ultra-cold atoms in optical lattices exemplifies the utility 
of these systems in studying the idealized lattice mod- 
els so central to condensed matter physics [IJ^. By in- 
creasing the depth of the lattice potential, an initially 
Bose-condensed system can undergo a transition from 
superfluid (SF) to Mott insulator (MI); careful experi- 
ments have pinpointed the critical lattice depth for this 
transition in 2D ^ and 3D [6,. The traditional Bose- 
Hubbard model describes homogeneous systems, but 
trapped ultra-cold gases are globally inhomogeneous, po- 
tentially containing multiple, spatially separated phases. 
For sufficiently large systems this inhomogeneity can be 
understood using the local density approximation (LDA) , 
where each region of the system is treated as being locally 
homogeneous. Rigol et al. [7] introduced a "universal 
state diagram" describing the full configuration of coex- 
isting, spatially separated SF and MI phases in harmoni- 
cally trapped systems. Here we present measurements of 
2D trapped systems and identify the transition from SF 
to MI as a function of lattice depth and atom number; 
the resulting experimental state diagram. Fig. [T] is in 
good agreement with the quantum Monte Carlo (QMC) 
predictions of Ref. [7], which go beyond the LDA. 

The BH hamiltonian models the physics of bosons in a 
lattice potential, here realized with ultra-cold *^Rb atoms 
in a 3D optical lattice. The homogeneous BH model 
includes only pair-wise on-site interactions and nearest- 
neighbor tunneling, parameterized by an interaction en- 
ergy U and a tunneling matrix element t. At zero tem- 
perature this model predicts the existence of a SF phase 
and MI phases with integer occupation n — 1,2,3... per 
lattice site, determined by the ratio U /t and the chem- 
ical potential /i. The importance of interactions is de- 
termined by U jt while the density is largely controlled 
by [Lit. For weak interactions (small U jt) the system 
is SF, while for U jt larger than a critical value {U/t)c 
the system can enter a MI phase. For U/t ^ {U/t)c, 
the phases alternate between SF and MI, increasing in 
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FIG. 1. Measured state diagram for a 2D Bose gas in the 
presence of harmonic confinement. The experimental transi- 
tion boundary (red ovals) was measured from f{U /t) traces 
like that of Fig. |4] for different N2D. Here, the size of the 
ovals denote the uncertainties in the measurement of p and 
U/t where we identify the phase transition. The data dis- 
played in Fig. |4] was taken along the green dashed line. The 
continuous blue line indicates the QMC prediction [J for the 
first appearance of a MI in the universal state diagram. The 
points are colored according to the side of the transition on 
which they are. 



density as /i increases [5]. 

The homogeneous BH model is not applicable to cur- 
rent trapped, ultra-cold atom experiments, owing to their 
harmonic trapping potential. The BH Hamiltonian for a 
lattice with period d superimposed on a symmetric har- 
monic trap is Ij 
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where b, is the creation operator of a boson at site i and 
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{i,j) constrains the sum to nearest neighbor tunnehng. 
The parameter e = muj'^dP /2 describes the harmonic po- 
tential where m is the atomic mass and v = 27ra; is the 
trap frequency. 

In the LDA, the third term of the Hamiltonian is as- 
sumed to be constant over an extended region producing 
a local chemical potential /i^ = fj, — ei^; under this ap- 
proximation the system's properties are computed using 
the homogeneous BH Hamiltonian. The LDA intuitively 
explains the evolution of a SF system into a nested col- 
lection of alternating SF and MI shells asU/t increases. 
Such structures were first observed using a magnetic res- 
onance imaging (MRI) approach for a 3D system [S], by 
measuring coUisional shifts 110], and more recently by 
direct imaging [11 . 

The inhomogeneity introduced by the trap leads, for 
sufficiently large e, to a breakdown of the LDA, both 
quantitatively and qualitatively. For example, the trap 
potential can increase the critical value {U/t)c at which 
a Mott insulator first appears [7]. This suggests that 
a proximity-like effect [12] stabilizes the SF component 
where MI was expected. 

For a, T — trapped system, the three parameters: 
U /t, e and the atom number N, fully specify the quan- 
tum state of the system, even with coexisting regions of 
SF and MI. Somewhat surprisingly, only two indepen- 
dent variables are sufficient [7]: U/t and a characteristic 
density p = Ne/t. By monitoring the dependence of con- 
densate fraction f on p and U/t [5], we experimentally 
measure the state diagram for 2D systems in an optical 
lattice (Fig.[l]). 

In our experiment we partition a 3D Bose-Einstein 
condensate (BEC) into an ensemble of nearly indepen- 
dent 2D systems using a ID optical lattice along z (see 
Fig. [2]). Additional optical lattices along x and y pro- 
vide the periodic landscape for the BH Hamiltonian in 
each 2D system. In this configuration a critical value 
{U/t)c = 15.8(20) [I3], for which the n = 1 MI first ap- 
pears in the ensemble system, was measured in Ref. 0, 
at higher temperatures than discussed here. That ensem- 
ble measurement, however, could not distinguish between 
2D systems with different p. To overcome the ensemble 
averaging, we developed an improved MRI approach to 
slice out a small subset of nearly identical 2D systems 
and measure their momentum distribution. In addition, 
we use matter wave focusing [14J to more accurately iden- 
tify the condensate, thereby reducing the measurement 
uncertainty in /. 

We prepare a = 2x 10^ atom ^^Rb BEC |T3], with no 
discernable thermal component, in the |F = l,mp = 1) 
state in a 3D harmonic trap with measured trapping 
frequencies = {23.2(5), 27.4(3), 42.8(9)} Hz. 

The trap arises from a combination of optical, magnetic 
and gravitational potentials. The BEC is confined at 
the intersection of a pair of 1064 nm laser beams, prop- 
agating along X and y, with waists (exp(— 2) radii) of 
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FIG. 2. a) Our 3D BEC is divided into 2D systems by a 
deep optical lattice in the presence of a linear magnetic field 
gradient, both aligned along the z direction. We use a MRI 
technique to selectively address a small subset of 2D systems, 
b) Matter-wave focusing is used to better resolve the SF phase 
of the 2D Bose gas; here we show a schematic of the focusing 
of a single 2D system after free evolution during TOF. c-d) We 
compare the measured atom distribution (approximating the 
momentum distribution) of the ensemble of « 60 2D systems 
without focusing (c) with that of the addressed 2D systems 
in the mp = sub-level with focusing (d), both with an x-y 
lattice at 9.5 Er. 



about 55 pm. The BEC is 620 pm above the zero of 
a quadrupole magnetic field. At the center of the BEC 
the magnetic field is Bq = 193 pT, corresponding to a 
Zeeman shift gp^Bo/h — 1.35 MHz. The magnetic po- 
tential, nearly linear along z with a gradient of 2.180(4) 
kHz / /im, almost exactly cancels gravity and adds a har- 
monic anti-trapping potential in the x-y plane for our 
|F = 1, mp = 1) atoms. 

We load the BEC into a 3D optical lattice at the inter- 
section of three pairs of linearly polarized nearly counter 
propagating laser beams from a A = 810 nm TiiSapphire 
laser [TB] . These beams form independent ID optical lat- 
tices along x, y and z. The z lattice is always set to a 
final depth of 24 En and partitions the 3D BEC into a set 
of « 60 2D systems; the depth of the x-y lattice ranges 
from to 20 En and determines the parameter U/t. To- 
gether all confining potentials determine e. The recoil en- 
ergy and momentum are En — h^kj^/2m = h x 3.4 kHz 
and kn — 2tt/X. The lattices are turned on from zero 
intensity in 100 ms with a half-Gaussian intensity ramp 
(rms width of 37 ms). This time scale was chosen to 
be adiabatic with respect to interactions and all relevant 
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single particle energy scales [51 [T7] . We measure lattice 
depth to within w 2% by pulsing on each lattice sep- 
arately for 4 to 6 /iS and observing the resulting atom 
diffraction [TH],[TS]. 

We implemented a MRI approach to selectively ad- 
dress a localized group of nearly identical 2D systems, as 
schematically illustrated in Fig. [2^. A radio-frequency 
(rf) magnetic field Brf, linearly polarized along x + ij 
transfers atoms from tof = 1 to mp — and mp = — 1- 
We choose B^f to maximize the transfer into mp — us- 
ing a 400 fjis Blackman pulse (perfect transfer to mp = 
is impossible for our 3 level system). The 2 kHz rms 
spectral width of this pulse, combined with the magnetic 
field gradient along z gives a 0.9 ^.m. rms spatial resolu- 
tion (« 2 lattice-sites), resulting in the extraction of 2D 
systems with up to 4000 atoms. 

Following the rf pulse, the lattice potentials are ramped 
off with exponentially decreasing ramps (400 /xs time con- 
stant) - nearly adiabatic with respect to single particle 
energy scales in the optical lattice - approximately map- 
ping the occupied crystal momentum states in the low- 
est Brillouin zone to free momentum states [20]. At the 
same time, we remove the optical dipole trap in < 1 /xs, 
and the magnetic field gradient along z in 3 ms; the 
atoms then expand for a 18.1 ms time-of-flight (TOF). 
During part of the TOF, a magnetic field gradient ap- 
proximately along y quickly separates the three mp com- 
ponents. We then detect the final spatial distribution of 
all three components using resonant absorption imaging, 
which gives the approximate momentum distribution of 
each spin component separately. The mp = distribu- 
tion directly measures the momentum composition of the 
nearly identical 2D systems selected by the rf pulse, vir- 
tually eliminating the inhomogeneous averaging that is 
present in the mp = 1 distribution (see Fig. [2j:). 

Absorption images after TOF can differ from the in 
situ momentum distributions for two primary reasons: a) 
interactions during TOF and b) finite TOF, here 18.1 ms. 
We mitigate each of these effects as follows: a) The al- 
ready weak interactions during TOF for the small num- 
ber of atoms transferred into the mp = state are fur- 
ther reduced by the rapid expansion along z after release 
from the tightly confining vertical lattice, b) We used a 
matter-wave focusing technique - a temporal atom lens 
- that "images" the in situ momentum distribution for a 
finite TOF [14] . To focus the atoms we increase the har- 
monic trapping frequency by a factor of about 3, by lin- 
early ramping the intensity of our 1064 nm dipole trap in 
200 /iS, and then holding for 400 /iS (during the rf pulse) 
just before TOF. 

Our 3D BEG has a 56.9(4) ^m Thomas-Fermi radius 
after TOF. When partitioned into an ensemble of 2D 
systems the radius decreases to 47.2(5) ^m. For the ex- 
tracted 2D systems the radius is 19.9(2) /im. Finally, fig- 
ure[2](c-d) illustrate the dramatic reduction to 10.5(2) /im 
when both interactions and finite size effects are mini- 
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FIG. 3. Density profile n{z) for; a) a 3D BEC and b) 
an ensemble of 2D systems. The atom number calculated 
from the in situ Thomas- Fermi radius — 8.2(2) /^m is 
TVtf ~ 1-8(4) X 10^. The vertical dashed lines indicate the 
Thomas-Fermi radius from our fit. Continuous lines show a 
fit to the tn situ ID density profile n(z). The temperature 
of the selected 2D systems (squares in b)) is displayed on 
the right axis, as a function of position along z. On average 
r= 15(3) nK. 



mized. 

We carefully calibrated atom number by measuring 
the in situ ID density profile n{z), of a 3D BEC using 
the MRI technique (see Figj3^, circles). The Thomas- 
Fermi radius — 8.2(2) fim gives atom number iVxp = 
1.8(4) X 10^; direct integration of n{z) gives Ni^t ~ 
1.89(5) X 10^; measurement of absorption by all atoms af- 
ter TOF gives iVabs = 1.90(5) x 10^. These measurements 
show that the combination of shot-to-shot number fluc- 
tuations and number measurement uncertainty is ^3%. 
We confirm this by loading the BEC into the ID opti- 
cal lattice along z, and again measuring n{z). We find 
that the density profile expands along z (Fig. |3]d, circles) 
but the integrated atom number A'int = 1.84(5) x 10^ 
remains constant. Figure |3j3 also shows the measured 
temperature T in a ID optical lattice as a function of z 
(squares). T = 15(3) nK is nearly uniform over all sig- 
nificantly occupied lattice sites, indicating that the 2D 
systems taken together are effectively in thermal equilib- 
rium. Loaded into a shallow lattice this corresponds to 
a temperature k^T = 0.9(2)t. 

In our experiment we set U jt using the x-y lattice 
depth and p by rf-selecting 2D systems with the desired 
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atom number from among the « 60 available systems. 
As a result, each measured momentum distribution cor- 
responds to a single point on the U/t — p plane and we 
use the condensate fraction / to distinguish between the 
SF and MI phases. 

We experimentally define / as the fraction of atoms 
in the sharp, focused feature in the momentum dis- 
tribution. To remove the broad background, present 
due to thermal effects and quantum depletion, including 
atoms in the Mott phase, we fit to the thermal distri- 
bution of non-interacting classical particles in a 2D sinu- 
soidal band [5]. We smooth the fit function in a region 
within 0.1 k^i of the edge of the Brillouin zone to account 
for non-adiabaticities in the lattice turn off near the band 
edge PU]. We exclude a disk with 0.16 kn radius around 
the condensate feature from the fit and identify the con- 
densate as the atoms that remain within the disk after 
subtracting the fit. 

Figure [4] shows the condensate fraction / as a func- 
tion of U/t for 2D systems with N2D ~ 3500, / up 
to 0.8, and an initial temperature T = 0.9(2)t, a factor 
of 2 lower than that reported in Ref. [5] where / < 0.4 
and T « 2t. f(U/t) exhibits three regions. For small 
U/t, f decreases (fit to a line) until at / « 0.12 and 
{U/t)c — 21(2) the slope changes markedly. We asso- 
ciate this feature with the first appearance of a MI and 
the subsequent decay (fit to a parabola) with the spatial 
growth of the MI regions. For U/t > 60 the condensate 
fraction is indistinguishable from zero. The critical point 
for appearance of MI (U/t)c — 21(2) is consistent with 
trapped system QMC {iU/t)c = 20.5 at p = 53) [7]. 
A similar analysis at p « 20 shows the first appear- 
ance of MI at {U/t)c = 19(2), consistent with past mea- 
surements {{U/t)c — 15.8(20)) [S], homogeneous system 
QMC calculations {{U/t)^ = 16.5) [HHH], and trapped 
system QMC {{U/t)^ = 20 at p = 20) [HlMl- 

Our measurements are summarized in Fig.[T] We sam- 
pled about 1300 images with p up to 100 and U/t up 
to 100 and for each image extracted the condensate 
fraction /. The points are colored according to the side 
of the transition on which they are, light grey symbols 
denote the SF region and dark grey symbols show points 
with some MI. Red ovals mark the experimental phase 
transition boundary, their widths denote the uncertain- 
ties in the measurement of p and U/t where we identify 
the phase transition. 

Figure [l] corresponds with the predicted state dia- 
gram [7]; the continuous curve shows the expected first 
appearance of a MI. The shift of this curve to larger U/t 
for increasing p, reproduced by the data, goes beyond the 
LDA prediction where all phase transition lines are verti- 
cal, independent of p, once MI shells form. The discrep- 
ancy for p < 15 is expected due to increased sensitivity 
to thermal effects at low density where the SF transition 
temperature is extremely low. 

During the preparation of this paper we learned of a 
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FIG. 4. Condensate fraction of a 2D Bose gas 
with ~ 3500 atoms measured through the SF to MI 

transition; red dashed curves are the fits described in the 
text. Insets (a-c) display the average momentum distribu- 
tion n{k) = [nx{k) ny{k)]/2 at different U/t, where Uxik) 
is the momentum distribution integrated over y and likewise 
for ny{k). We distinguish the formation of the first and sec- 
ond Mott regions at U/t = 21(2) and U/t ~ 55 respectively. 
For U/t > 60 the condensate fraction is indistinguishable from 
zero. 
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